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We investigate theoretically the quantum oscillator-like states recently observed experimentally 
in polariton condensates (Nat. Phys. 8, 190 (2012)). We consider a complex Gross-Pitaevskii type 
model which includes the effects of self-interactions, and creation and decay of exciton-polaritons. 

We develop a perturbation theory for approximate solutions to this non-equilibrium condensate 
model and compare the results with numerically calculated solutions for both repulsive and attractive 
polariton-polariton interactions. While the nonlinearity has a weak effect on the mode selection their 
density profiles are modified at moderate gain strengths and becomes more dominant when a very 
large gain of polaritons implies an extended cloud with high condensate densities. Finally, we 
identify the relation of the observed patterns to the input pump configuration, and suggest this may 
serve as a generalized NOR gate in the tradition of optical computing. 

PACS numbers: 03.65.-w, 02.60.Cb, 03.65.Sq, 03.75.-b 


INTRODUCTION 

Bose-Einstein condensates (BECs) are prime examples 
of nonlinear quantum many-body systems, because of in¬ 
teractions between coherent quantum particles. They are 
regarded as a macroscopic quantum phenomenon due to 
their very high coherence and long-range order of all par¬ 
ticles constituting this state of matter, a quantum state 
that typically is formed at very low temperatures m- 
Many setups of atomic BEC are for fixed particle num¬ 
bers and can carry various topologically stable excita¬ 
tions ranging from solitons to giant vortices 0 ehei. 
More recently efforts have been made towards the real¬ 
ization of coherent non-equilibrium many-body systems 
that exchange particles and/or energy while at the same 
time obeying BEC properties e.g. in solid state systems 
on ns. One such system is the polariton quasiparti¬ 
cle Bose-Einstein condensate mummm, which can 
be effectively described by a non-conservative nonlinear 
Schrodinger-type equation El US El formally generaliz¬ 
ing the so called Gross-Pitaevskii equation (GPE) - the 
model for the condensate wave function of dilute and 
weakly interacting atomic BEC [imim. 

Bose-condensed polaritons have become feasible in 
semiconductor microcavities in the last decade mMm 
m- Polaritons in semiconductor microcavities are su¬ 
perpositions between excitons (electron-hole pairs) and 
cavity photons, and as such the combination of the com¬ 
ponents determines the particle statistics [2j [15] . These 
quasiparticles in a dilute regime show distinctive prop¬ 
erties such as Bose-Einstein condensation, superfluid¬ 
ity, a finite nonzero speed of sound due to the non¬ 
linearity and the emergence of elementary topologically 
stable excitations such as solitons or quantum vortices 
[2l EUm ng HZ] . In their very nature polariton BECs are 
non-equilibrium systems - polaritons can be created via 


a local light field and decay after 1-100 ps mainly due to 
leaking of the cavity photons mmm- The decay rate 
significantly depends on the quality of confinement of the 
light field [18] within the microcavity and in this work it 
is assumed to be constant over space. From a techno¬ 
logical and experimental point of view polariton BECs 
have several practical advantages over atomic BECs. For 
instance, BECs of dilute atoms must be realized at tem¬ 
peratures within the nano-Kelvin range [MS], whereas 
polariton condensates can be observed even at room tem¬ 
perature in polymers [22] and in organic materials [23] . 
This is a consequence of the fact that polaritons have a 
very light mass due to their inherent coupling to photons 
0, and so have a much higher critical Bose condensa¬ 
tion temperature. Restrictions on the polariton density 
however have to be taken into account to ensure their 
bosonic character is maintained [3]. 

In many settings a polariton condensate is an effec¬ 
tively 2D coherent many-body system that is tightly con¬ 
fined with regards to the optical axis of the semiconduc¬ 
tor microcavity [HEEIJEIE]. Depending on additional 
constraints the condensate may become essentially ID 
mm- The elementary excitations within a ID system, 
with repulsive polariton-polariton interactions, are dark 
and grey solitons mm- Grey solitons are stable excita¬ 
tions which move over long distances without changing 
their form and in principle with elastic scattering prop¬ 
erties, but due to the non-equilibrium character obey a 
finite lifetime [25]. Dark solitons are the stationary spe¬ 
cial case of the grey soliton when the density depletes 
to zero at its minimum and with a tt phase shift at the 
nodal point [4 . 

In this work we are particularly interested in station¬ 
ary and stable patterns within a ID polariton conden¬ 
sate, that once excited are similar to dark soliton jj3| ar¬ 
rays [T6] with respect to singularities at localized points - 
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the oscillator modes of a harmonically trapped quantum 
state. Experimental observations of oscillator modes in 
polariton condensates have been made in [23123, but a 
thorough theoretical description is lacking. The exper¬ 
imentally observed non-equilibrium quantum oscillator 
states were formed between two spatially separated pump 
spots of incoherently pumped polaritons, corresponding 
to two localized laser beams shining on the semiconduc¬ 
tor microcavity and positioned within the 2D plane. Be¬ 
tween the pump spots stationary spatial patterns were 
formed similar to those observed for a quantum harmonic 
oscillator, hence suggesting an effective ID problem. 

The spatial patterns can be understood in a dynam¬ 
ical picture as follows. The polaritons generated at the 
spots move from the source peaks in all directions, be¬ 
cause interactions of pump and reservoir with polaritons 
are repulsive [24]. So in between the pump spots the 
bulk of polaritons are present and thus the condensate 
and its patterns form predominantly there. Depending 
on the pump strength, and the repulsiveness from the 
pump/reservoir, different modes can be observed; start¬ 
ing from the ground state of approximately gaussian form 
(when nonlinearity plays a minor role), excited states be¬ 
come feasible. Experimentally only the case for repulsive 
self-interactions has been considered so far. In this work 
we predict - based on the recent experimental effort [28 - 
the states for attractive self-interactions as well and their 
dependence on the gain. To elucidate the role of nonlin¬ 
ear self-interactions for the pattern formation in quasi- 
harmonically trapped polariton condensates we vary be¬ 
tween attractive and repulsive interactions and the non¬ 
linear character of polariton BEC will become dominant 
in a Thomas-Fermi (TF) type regime, i.e. when conden¬ 
sate densities are very high and pumping is strong. 

Our paper is organized as follows. First we will in¬ 
troduce the general model describing the polariton con¬ 
densate in 2D and a simplified model for the effective 
ID investigations. Then approximate analytical expres¬ 
sions for the condensate wave functions are presented and 
compared with numerical simulations. Subsequently the 
dependence of the eigenvalues on pumping strength and 
the nonlinearity is presented. Finally we introduce a gen¬ 
eralized optical NOR device and summarize our findings. 


POLARITON CONDENSATE MODEL 

We describe the order parameter of the polariton 
condensate by the state equation in the spin coherent case 
m is, which has proven to be an accurate description 


of recent experiments pa is], 

ia <*=\-^ + v< N + 

+ gpP + ai\i>\ 2 + ^(gpN (1) 

where m* is the effective mass of the condensed polari¬ 
tons and P describes the geometry and magnitude of the 
pump spots, which are induced by linearly polarized laser 
beams [ 12117 ]. The polariton reservoir is approximately 
given by mm 



In Eq. [l] ye denotes the globally constant decay rate 
of condensed polaritons, 7 r the relaxation of the reser¬ 
voir, (3 is the scattering rate of polaritons into the con¬ 
densate, aq denotes the polariton-polariton interactions 
within the condensate, gp is the repulsive interaction of 
the condensate with the pump and gp the repulsive inter¬ 
actions of the condensate with the reservoir. Parameter 
values are chosen in accordance to most recent experi¬ 
ments nansi, i.e. effective mass m* = 5x10 5 mo where 
mo is the mass of the free electron, aq = 0.001*/im 2 , 
9r = 0.022 • ps~ x • pm 2 , gp = 0.07 • ps -1 • pm 2 , (3 = 
0.05 • ps~ x • pm 2 , 7 r = 10 - ps -1 , Rr = 0.06 - ps -1 • pm 2 , 
7 c = 0.556 • ps _ 1 . The pump geometry is assumed to be 

P(x,y) = A^exp (-0.073 • ((x + aq) 2 + y 2 )) + 

T exp (-0.073 • ((x - x 2 ) 2 + y 2 )) ^ (3) 

with A being the amplitude and the position of the two 
gaussian pump spots is at x 2 and aq respectively, which 
models the gain due to local illumination of the semicon¬ 
ductor microcavity. If we neglect spatial extension in y 
we simply write P(x). We note that the exact mathemat¬ 
ical form does not significantly change the physics under 
consideration and is an approximation of the actual spot 
form used in recent experiments Hal m]. 

Taking this into account, and as outlined in the next 
section in more detail, we can effectively reduce 0 to a 
simplified model for the order parameter i/j of the polari¬ 
ton condensate given by [30 

= ( K ~ n ^+v(r)+uw\ 2 +i( 7eff-r|V’| 2 ))v>- 

( 4 ) 

Here U denotes the self-interaction strength, 7 e ff the lin¬ 
ear gain and T the nonlinear loss. The repulsive effect 
on the polaritons of the pump P and the reservoir A, 
defined in 0, is assumed to be approximated by V. For 




3 


X 



FIG. 1. Effective harmonic confining potential (dashed line) 
with k = 0.05 vs. the pump profile P (dotted line) § for 
A = 5, X 2 = 10 iim and x\ — —10 [im. A harmonic oscillator 
mode forms within the trapping (continuous line). 

our analysis we are interested in stationary states and so 
set ihdt'ip = H'lp- As a simplification, we neglect the trans¬ 
verse spatial dimension, and so write ip(x, y) = 'ip(x), and 
rescale @ as in [30], i.e. V(x) = (Jiw/2)( k x 2 /l 2 ) (however 
without effective transverse confinement), where uj is the 
oscillator frequency and l = %Jh/ (m*uj) the correspond¬ 
ing length. Expressing lengths in units of l and energy in 
units of Hu, and by substituting ip y/ftLVjJZU)ip, the 
result is 

(dxx - X 2 + (<Ti +ia 2 )\ip\ 2 - inPjip = -MiV’ (5) 

with ni = (|£), M 2 = 27 e ff/ (hu>), <j\ = ±1 (+ corre- 
sponds to attractive and — to repulsive self-interactions), 
(72 = T/U and V{x) = x 2 is the harmonic potential 
due to the approximation of the pumping spots left and 
right of the condensate and the repulsive reservoir. We 
note that is a nonlinear eigenvalue problem with real¬ 
valued eigenvalue /ii, which we will determine numeri¬ 
cally and analytically. 

Let us identify the Carusotto-Wouters model as ap¬ 
proximated in [H|2l|29] with the Keeling-Berloff model 
m by identifying the effective potential terms by 

guN + gpP + a^l 2 ^V(r) + U\^\ 2 (6) 

with V = kx 2 and the effective gain and loss terms 

( 7 i?V- 7c )/ 2 o ( 7eff -rH 2 ) (7) 

with N = ^(1 — y 1| 2 ). To approximately match the 
potential terms, we neglect the density dependent part 
of N (as it does not significantly affect the observed pat¬ 
terns) and we use a harmonic function that passes the 
maxima of the two Gaussian spots of |3|. In Fig. [l] we 
show the correspondence between the pump spots and 
the harmonic approximation and remark that k ~ A. To 
fit parameters in |6| we note that the density dependent 
potential term can simply be matched aq = U. 

To identify the pump terms we set P — 7 c = 2 y e ff, 
where P = l/(\x% — X 2 I) f** Pdx = 0.33 • A is the 


average value of gain between both peaks assumed to 
be at x\ = —10, x 2 = 10 and 2T = P/^r. Con¬ 
sequently the density dependent pump parameter be¬ 
comes (J 2 = T/U = 2.5 while the linear parameter is 
/i 2 = (0.165 • A - 0.278)/fc with k = m*uj 2 /2 = 0.092 • A 
as in m fixing the potential strength uj via the choice 
of A. This identification yields the simplified equation 
([ 5 ]) once the corresponding rescaling has been applied. 
Recall that sgn(cri) depends on the sign of the scattering 
length between condensed polaritons. 

By tuning the properties of the pump (in particular the 
geometry and pumping strength), and also having con¬ 
trol over the nonlinear interaction strength oq through 
the use of Feshbach resonances [28] , we see that the gain 
and loss coefficients ((72 and ^ 2 ) can be varied widely. 
In the remainder of this work we will set the density 
dependent loss rate in rescaled units (72 = 0.3 and con¬ 
sider values for the pumping strength /ii from fi^ — 0 to 
li 2 = 1, values that are in accordance with Wi- We re¬ 
mark that generally loss rates for polariton condensates 
depend significantly on the actual experimental setup, 
partly because lifetimes vary from lps to lOOps and cor¬ 
respond to the quality of confinement within the micro¬ 
cavity 0 mia E3 ng, so even a quasi-equilibrium scheme 
is feasible in polariton condensates m ■ This allows us 
to examine the equilibrium limit in our nonlinear sys¬ 
tem seamlessly (/12 = 0), and explore the nature of the 
system when gain and loss are of similar magnitudes but 
less than or equal to the magnitude of the nonlinear term. 
Setting cr 1 — zbl means an increase of real-valued ‘non¬ 
linear’ behavior than if we used the initially introduced 
‘natural’ parameters, which is justified in physical terms 
by utilizing Feshbach resonances mmm- 

THE POLARITON HARMONIC OSCILLATOR 

First we present results for simulations of the unre¬ 
duced state equation Eq. [l] The pump geometry is 
assumed to be (|3|. In Fig. [2] (a) we present the den¬ 
sity profile IV’I 2 showing the harmonic oscillator ground 
state in-between two pump spots for A = 8.5 roughly 
approximated by a gaussian. Fig. [ 2 ] (b) shows the first 
excited harmonic oscillator-type state for A = 11.5 obey¬ 
ing two distinct peaks of condensate density. Finally Fig. 
[ 2 ] (c) shows the second excited harmonic oscillator-type 
state for A = 35 obeying three distinct peaks of conden¬ 
sate density. Increasing the amplitude allows us to access 
higher order excited states. 

Although we mainly aim to explain the patterns ob¬ 
served in [26], we note that the polariton harmonic oscil¬ 
lator can itself be implemented within a long ID nanowire 
similar to PH CHS i-e. a microcavity with very strong 
transverse confinement. Additionally an external in¬ 
plane trap can be imposed upon the effectively ID con¬ 
densate as it has been experimentally shown for a sample 
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FIG. 2. Density profile |'0(x,y)| 2 of the harmonic oscillator 
ground state in a 2D setting for an area of 10 x 10 jim is shown 
in (a). In (b) \ip(x,y) | 2 of the first excited harmonic oscillator 
state in a 2D setting for an area of 10 X 10 jim and in (c) the 
density profile of the second excited harmonic oscillator state 
in a 2D setting for an area of 10 x 10 \im. m 


consisting of three sets of four GaAs/AlAs quantum wells 
embedded in a GaAs/AlGaAs microcavity in [35, 36j. 
On the other hand for the pump laser beams there is in 
principle no restriction to the size of the area they are 
illuminating m and thus they can be assumed to be 
extended over a wide range, including the area were the 
condensate forms within the harmonic trap. So an almost 
homogeneous pump distribution can be provided on the 
ID geometry on top of a repulsive and approximately 
harmonic trap as described by the simplified model <©• 


PERTURBATION THEORY 


We now derive analytical expressions approximately 
solving (§. Considering <§ in the parameter limit 
ai , <J 2 , /i 2 0 we get the well-known quantum harmonic 


oscillator equation 


(pxx % ) — E n (/) n , 

(8) 

which has Hermite-Gauss polynomials as 

eigensolution 

fl b n (x) = c n e- x2/2 H n (x ) 

(9) 

with c n = (2 n n!y / 7r) _1 / 2 , where 


2 d n 2 



(10) 


with corresponding eigenvalues E n = 1 + 2n with n G 
{0,1,2,3,...}. These Her mite-Gauss polynomials are an 
orthonormal basis (ONB) of a Hilbert space. On the 
other hand any wave function in a Hilbert space has an 
expansion in terms of Hermite polynomials of the form 

oo 

tp( x ) = 'E Bn< ^ n ^’ 

n =0 

with complex coefficients and real-valued functions </> n (x) 
providing the ONB. By inserting this expansion in (J5| 
and using we get 


oo / OO OO \ 

y ^ ( E n + (<Ti + i(J2) 'y ^ B k cj) k y ^ Bi (f)i j B n cf) n = 

n =0 \ k=0 1=0 J 

oo 

= ( — /H + ^ 2 ) B n cj) n . (12) 

n =0 

By multiplying from the l.h.s. with </> m and integrating 
over the whole space we obtain the formula 


(^"l H - io' 2 ) y ^ BfcBi Bj^yn k i n — 
k,l,n=0 

= (Em fi\ -\- B rn ^ (13) 

where we have introduced the abbreviation V m kin — 

J> rn (\) k (\)i(\) n dx. 

We assume that for the j th nonlinear mode the largest 
coefficient in our expansion © is Bj. We take B n< j = 0 
and assume that only basis wavefunctions of the same 
symmetry of the nonlinear mode contribute in the expan¬ 
sion, thus we take R^fc+i = 0 (where k G N). Therefore 
at leading order ( fTdj ) becomes: 

(a 1 E icr 2 )\ \Bj\ 2 Vjjjj = (Ej — /ii + i/i 2 ). (14) 


Solving simultaneously the real and imaginary parts of 

(15) 


(14) gives us an expression for the eigenvalue /ii , 

01 M 2 


/i\ Ej 




and for \B« 


\Bj \ 2 = —rr 

a 2 Vj 


(16) 


3,3,3,3 
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In our comparison with the numerical results in the next 
section we will use only these leading order results as 
they are relatively simple and easy to calculate, and im¬ 
mediately give us an indication of how “nonlinear” the 
rescaled wavefunction is. The eigenvalue (15) depends 
linearly on the nonlinear coefficient, but the amplitude 
of the rescaled wavefunction, given by the square root 
of (16), depends just on the gain and loss parameters, 
in a way which is in accordance with the results due 
to a non-conservative Lagrangian formalism presented in 
m o, when back scaling of the perturbation series to 
the wave function satisfying Q is made. Thus if we see 
that our numerical solution agrees well with this leading 
order term then we can conclude that the nonlinearity 
plays only a minor role for the rescaled condensate wave 
function in this particular parameter regime. 

We see that at leading order we have no information 
about the phase of the wavefunction. To proceed fur¬ 
ther and find simple expressions for the higher order 
contributions we assume that Bj is purely real and that 
B;, Bj-^-2k' This allows us to find the following ex¬ 

pression for the higher-order complex coefficients: 


B 


'j+2k 


(Ji + icr 2 


Ej+2k ~ /M + 


B jV 2k 




(17) 


Note that despite considering a non-equilibrium system 
the condensate wave function satisfies the normalization 
condition f|^| 2 dx = \B^\ 2 = c < oc and c cor¬ 

responds to the number of condensate particles similarly 
to the case of fixed particle numbers [4]. 

Thomas-Fermi type regime. - Starting from (J5| we con¬ 
sider the regime where the kinetic energy becomes neg¬ 
ligible compared to the external potential, nonlinearities 
and linear parameters, i.e. 


-x 2 + (<Ti + ia 2 )\ip \ 2 - in 2 )v> - -MiV’ 


(18) 


which to the leading order has the formal solution 


r b (x) = ±< 


~ni + in 2 + x 1 

(Xi + io 2 


(19) 


Expression (19) is however not integrable for arbitrary 


parameters. Besides we refer to [SHU for a rigorous dis¬ 
cussion of the Thomas-Fermi approximation in conserva¬ 
tive GP theory. In contrast to the previous perturbation 


theory, where the coefficients (16) and thus the conden¬ 


sate wave functions dependent solely on the parameters 
ji 2 and (72, the TF regime significantly depends on the 


nonlinearity (j\ and the energy fii as formula (19) shows. 


J 


The validity of the leading-order analytical prediction 


COMPARISON BETWEEN NUMERICAL AND 
ANALYTICAL RESULTS 


We now proceed to test our leading order perturba¬ 
tion theory by comparing it with numerically calculated 
results. We use a modified squared-operator method [38] 
to find numerical solutions to the equation <§• 

In Figure [3] we show the first three nonlinear modes for 
both repulsive (Fig. [3](a-c)) and attractive (Fig. §d-e)) 
interparticle interactions for the case ii 2 = 0.2. The nu¬ 
merical solutions are given by the thick lines, with real 
and imaginary parts corresponding to the dashed and 
dash-dotted lines respectively, while the density is given 
by the solid line. We see that, as in the case of the single¬ 
component atomic BEC m, the jth mode has j—1 nodes 
of zero density where a tt phase change occurs. However, 
unlike the atomic BEC case, the polariton condensate has 
a spatially dependent phase away from the phase singu¬ 
larity (evident through the spatially varying ratio of real 
to imaginary part in the wavefunction). This variable 
phase indicates that particle flow must occur in the con¬ 
densate, with particles moving from regions where gain 
dominates, to regions where loss dominates. 


We can compare the numerical results directly with the 
analytical expansion (11) assuming only the leading or¬ 
der term is present, as given by Eqs. (15) and (16). The 
analytical results are given by the thin solid lines in Fig. 
[3] More precisely we can only compare the analytically 
predicted density with the numerical result, however we 
include also the real part and imaginary part assuming 
that the imaginary part is zero. At this leading order 
the nonlinearity plays no role in the form of the ana¬ 
lytical wavefunction, so it is identical for the attractive 
and repulsive cases. Nevertheless, we see that it agrees 
remarkably well with the numerical solutions. We see 
that in the repulsive case the true solution lies outside 
the analytical result, while the inverse is true for the at¬ 
tractive case, as expected. Interestingly we see that the 
agreement improves as we go to higher modes, suggest¬ 
ing the ground state is “more nonlinear” in the sense 
that it deviates more strongly from the associated linear 
wavefunction. 


In Figure [4] we compare the analytical and numerical 
results at higher gain, given by fi 2 — 0.5. We see that 
the pattern observed for fi 2 = 0.2 becomes more evident, 
with the ground state showing the strongest deviation 
from the linear prediction. The higher-order modes con¬ 
tinue to show remarkably good agreement. 


can be seen at a glance by comparing the predicted eigen- 
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FIG. 3. Comparison between numerical (thick lines) and analytical predictions (thin lines) for repulsive interactions (a-c) and 
attractive interactions (d-e) at gain strength fi 2 = 0.2. (a,d) Ground state; (b,e) First excited state; (c,f) Second excited state. 
Numerical results follow format: density - thick solid lines; real part - thick dashed lines; imaginary part - thick dash-dotted 
lines. We see increasingly good agreement between the numerical and analytical predictions as we go to higher modes. 


values ( p~5] ) with those found numerically. We see in Fig¬ 
ure [5] the dependence of the eigenvalues on gain parame¬ 
ter fi 2 for the first three modes (pairs of lines are arranged 
upper to lower and correspond to second, first zeroth 
modes respectively), for both repulsive (solid) and at¬ 
tractive (dashed) interactions. The analytical predictions 
are given by the associated thin lines. We see increas¬ 
ingly good agreement with mode order, with the second 
mode showing close agreement even up to (i^ = 1. The 
ground state on the other hand is showing a discernible 
discrepancy by fi 2 = 0.5. The vertical dash-dotted lines 
correspond to the values of /i 2 chosen to show the sta¬ 
tionary states in Figs. [3]and[4j Overall we see that the 
features of the stationary states are dictated predomi¬ 
nantly by the gain and loss coefficients, rather than the 
interparticle interactions. 

NOR GATE AND A nNOR GATE 

The pattern formation within a ID nanowire can have 
the following technological application as NOR gate. A 
NOR gate [39] is defined by the logic that two or more 
zero inputs yield an output 1 in a boolean logic system. 
Here we can identify each pump spot © with 0 or 1 
by defining that, if the laser beam is turned on, i / 0, 
the input corresponds to 0 and if not it corresponds to 
1. The output will be the quantum harmonic oscillator 


patterns observed, iff both pump spots are equal 0. So we 
can identify all the harmonic oscillator patterns with an 
output 1, and thus we have established a simple feasible 
NOR gate by means of topologically stable excitations 
carried by a non-equilibrium polariton condensate within 
a semiconductor microcavity which stand in the tradition 
of optical computing P3- 

This NOR gate can be generalized in the sense that we 
do not solely distinguish between ”on” corresponding to 
0 and ”off” associated with 1 of both pump spots, but 
in addition include the experimental and numerical ob¬ 
servation that different pump strengths lead to different 
numbers of density lobes. Therefore we can enumerate 
the logical Os leading to different lobes states, i.e. ground 
state, first excited state, second excited state and so on. 
That is, if the pump spots have a certain amplitude A it 
leads to patterns of say B(A) : M —>> N lobes. However 
when only one pump spot is turned on, corresponding to 
a logical input, 0 and 1, or boths pumps are off, corre¬ 
sponding to 0 and 0 input, no lobes patterns form H3- 

CONCLUSIONS 

Motivated by a recent experiment we considered a 
mean field model for polariton condensates to reproduce 
the observed patterns. Starting from this we could iden¬ 
tify the general model with that of a simpler nonlin- 
































FIG. 4. Comparison between numerical (thick lines) and analytical predictions (thin lines) for repulsive interactions (a-c) and 
attractive interactions (d-e) at gain strength fi 2 = 0.5. (a,d) Ground state; (b,e) First excited state; (c,f) Second excited state. 
Format of the lines is the same as in Fig. [3] We see the ground state shows significant discrepancies with the analytical 
predictions, however for the higher modes the agreement is still very good. 
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FIG. 5. Family of eigenvalues, i.e. first, second and third 
order modes as function of pumping strength / 12 . Solid lines 
represent the numerically calculated eigenvalues for the repul¬ 
sive case and correspondingly the dotted lines are the analyt¬ 
ical prediction. For attractive condensates the dashed lines 
show numerically generated predictions and dotted lines are 
the analytical predictions. 


ear non-equilibrium quantum harmonic oscillator. For 
this effective theory we have presented a perturbative 
approach and so generalized the quantum harmonic os¬ 
cillator theory to include the non-conservative character 
of the non-equilibrium polariton condensate. We intro¬ 
duced the family of excited states starting from the gaus- 


sian ground state to the n node non-equilibrium quantum 
harmonic oscillator states. The analytical eigensolutions 
were compared to corresponding numerical simulations 
and showed excellent agreement. Due to the recent ex¬ 
perimental accessibility of Feshbach resonances in polari¬ 
ton systems we further predicted the scenario of attrac¬ 
tive self-interactions and outlined their role in nonlin¬ 
ear pattern formation. Particularly the modification of 
the density formation was explicitly illustrated. In the 
course of our analysis we showed that the energy eigen¬ 
values depend on the pumping strength fi 2 linearly to the 
leading order and conclude to provide a thorough expla¬ 
nation for the experimentally observed patterns. Finally 
we discussed the possibility for a generalized optical NOR 
device in polariton condensates utilizing the underlying 
pattern formation logic. 


ACKNOWLEDGEMENTS 

FP has been financially supported through his EPSRC 
doctoral prize fellowship at the University of Cambridge. 
FP likes to thank Natasha Berloff, Hamid Ohadi and 
Alexander Dreismann for stimulating discussions. 


* florian.pinsker@gmail.com 
t t.alexander@unsw.edu.au 

[1] A. Einstein, Sitzungsber. Preuss. Akad. Wiss., 3 (1925). 



























[2] H. Deng, H. Haug and Y. Yamamoto, Rev. Mod. Phys. 
82, 1489 (2010). 

[3] A. J. Leggett, Rev. Mod. Phys. 73, 307 (2001). 

[4] L.P. Pitaevskii and S. Stringari, Bose-Einstein Conden¬ 
sation, Clarendon, Oxford, (2003). 

[5] F. Pinsker, Excitations in superfluids of atoms and 
polaritons , Doctoral Thesis, University of Cambridge, 
(2014). 

[6] I. Carusotto and C. Ciuti, Rev. Mod. Phys. 85, 299 
(2013). 

[7] F. Pinsker and N.G. Berloff, Phys. Rev. A 89, 053605 
(2014); arXiv:1401.1517 (2014). 

[8] M. Correggi et al. , Phys. Rev. A 84, 053614 (2011), 
arXiv:1107.5160 (2011). 

[9] M. Correggi et al. , The European Physical Journal Spe¬ 
cial Topics 217, 183-188 (2013), arXiv:1210.7134 (2012). 

[10] F. Pinsker, N.G. Berloff and V.M. Perez-Garcfa, Phys. 
Rev. A 87, 053624 (2013); arXiv: 1305.4097 (2013). 

[11] J. Kasprzak et al. , Nature 443, 409-414 (2006). 

[12] A. Baas et al ., Phys. Rev. Lett. 96, 17640 (2006). 

[13] F. Pinsker, arXiv:1501.01951 (2015). 

[14] F. Pinsker and H. Flayac, arxiv:... (2015). 

[15] N. G. Berloff and J. Keeling, Physics of Quantum Flu¬ 
ids, 177, pp. 19-38, Springer Berlin Heidelberg, Germany 
(2013). 

[16] F. Pinsker and H. Flayac, Phys. Rev. Lett. 112, 140405 
(2014); arXiv: 1310.7500 (2013). 

[17] F. Pinsker, Int. J. Mod. Phys. C DOI: 

10.1142/S0129183115500825 (2015), arXiv: 1305.4088 

(2013). 

[18] R. Nelsen, et al. , arXiv: 1209.4573 (2012). 

[19] M. H. Anderson et al ., Science 269, 198 (1995). 

[20] K. B. Davis et al. , Phys. Rev. Lett. 75, 3969 (1995). 

[21] C. C. Bradley et al. , Phys. Rev. Lett. 75, 1687 (1995). 

[22] J. D. Plumhof, et al., Nature Materials, 

doi: 10.1038/nmat3825 (2013). 

[23] S. Kena-Cohen and S. R. Forrest, Nature Photonics 4, 
371-375 (2010). 


[24] A. Dreismann et al, Proceedings of the National 
academy of Sciences of the United States, 111 , No. 24, 
8770-8775 doi: 10.1073/pnas.l401988111 (2014). 

[25] Y. Xue and M. Matuszewski, Phys. Rev. Lett. 112 , 
216401 (2014). 

[26] G. Tosi, et al., Nat. Phys. 8, 190-194 (2012). 

[27] A. Askitopoulos, T. C. H. Liew, H. Ohadi, Z. Hatzopou- 
los, P. G. Savvidis, P. G. Lagoudakis, arXiv: 1411.4579 
(2014). 

[28] N. Takemura, S. Trebaol, M. Wouters, M. T. Portella- 
Oberli and B. Deveaud, Nature Physics 10, 500 - 504 
(2014). 

[29] M. Wouters and I. Carusotto, Phys. Rev. Lett. 99, 
140402 (2007). 

[30] J. Keeling and N.G. Berloff, Phys. Rev. Lett. 100, 250401 
(2008). 

[31] Y. S. Kivshar, T. J. Alexander, S. K. Turitsyn, Phys. 
Lett. A 278, pp. 225-230 (2001). 

[32] M. A. Karim, A. A. Awwal, Optical Computing: An In¬ 
troduction, John Wiley and Sons, Inc. New York, NY, 
USA (1992). 

[33] E. Kammann, Doctoral thesis, University of Southamp¬ 
ton (2013). 

[34] The time evolution has been induced via a 4th order 
Runge-Kutta method and a 4 order finite differences 
scheme. We used a computational window of 226 space 
units in x and along the transverse directions y. Typi¬ 
cal space steps are about Ax = 0.1 and the time steps 
At~1.5xl0 -4 . 

[35] R. Balili et al., Applied Physics Letters 88, 031110 
(2006). 

[36] B. Nelsen et al, Journal of Applied Physics 105, 122414 
(2009). 

[37] private discussion with Dr Hamid Ohadi, Cavendish Lab¬ 
oratory, UK 

[38] J. Yang and T.I. Lakoba, Stud. Applied Math. 118, 153 
(2007). 

[39] M.M. Mano, and C. R. Kime, Logic and Computer De¬ 
sign Fundamentals, Third Edition. Prentice Hall, p. 73 
(2004). 

[40] E. Kammann et al, Phys. Rev. Lett. 109, 036404 (2012). 


